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Abstract. It is generally accepted that the atmospheres of cool/lukewarm stars of spectral 
types A and later are described well by LTE model atmospheres, while the O-type stars require 
a detailed treatment of NLTE effects. Here model atmosphere structures, spectral energy 
distributions and synthetic spectra computed with Atlas9/Synthe and Tlusty/Synspec, 
and results from a hybrid method combining LTE atmospheres and NLTE line-formation 
with Detail/Surface are compared. Their ability to reproduce observations for effective 
temperatures between 15 000 and 35 000 K are verified. Strengths and weaknesses of the different 
approaches are identified. Recommendations are made as to how to improve the models in 
order to derive unbiased stellar parameters and chemical abundances in future applications, 
with special emphasis on Gaia science. 



1. Introduction 

Quantitative stellar spectroscopy employs model atmospheres that are based on many approxi- 
mations and simplifications. One of the fundamental decisions to be made for the model 
construction concerns the description of the thermodynamic state of the atmospheric plasma. 
Either local thermodynamic equilibrium (LTE) may be assumed, or deviations from it may be 
allowed for (NLTE), i.e. taking the mutual interaction of radiation and matter into account. 

It is generally accepted that the atmospheres of cool/lukewarm stars of spectral types A and 
later are described well by LTE model atmospheres, while O stars require a detailed treatment 
of NLTE effects due to the high energy density of their radiation field. Both approaches are 
being followed in the literature for analyses of B-type stars at present, with variable success. 

The present work tests the predictive power of common LTE and NLTE model atmosphere 
and line-formation codes in this transition region, concentrating on main-sequence stars of 
effective temperatures between 15 000 and 35 000 K (corresponding to spectral types of about B4 
to 09). Three approaches are investigated: pure LTE as implemented in the Atlas9/Synthe 
suite of codes [TJ [2], pure NLTE with Tlusty/Synspec [3J, and a hybrid approach based 
on LTE atmospheres (Atlas9) and subsequent NLTE line-formation computations with 
Detail/Surface [H We deviate from the usual 'model vs. model' comparison - which 
may be misleading as a priori it is not clear as to which model is more realistic - by testing 
observational constraints wherever possible. 
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Figure 1. Near-IR spectrum of HR 1861 (Bl IV), observed with Foces on the Calar Alto 2.2m 
telescope. The most important lines in the Gaia RVS spectral range (grey region) are identified. 

2. The Gaia perspective for investigations of OB-type stars 

The legacy of the Gaia mission in terms of OB-type stars will mainly lie in the determination 
of highly precise parallaxes, proper motions and broad-band space photometry. In contrast 
to the cooler stars, little information besides radial and projected rotational velocities will be 
derivable from the spectra obtained with the on-board radial- velocity spectrometer (RVS). The 
reason for this is the scarcity of spectral lines in the wavelength range of the RVS, which 
essentially comprises the Paschen series members P13 to Pi6 in this kind of star, see Fig. [TJ 
The determination of accurate and precise stellar parameters and elemental abundances will 
therefore have to rely on ground-based observations, e.g. within the Gaia-ESO survey [BJ. 



3. Models: assumptions & implementations 

Quantitative spectroscopy of early-type main-sequence stars and giants in the effective tem- 
perature range of interest, 15 000 < T c g < 35 000 K, can rely on classical model atmosphere^ 
Basic assumptions for the model construction are chemical homogeneity, stationarity, plane- 
parallel geometry and hydrostatic and radiative equilibrium, while the thermodynamic state of 
the atmospheric plasma can be determined in LTE or NLTE. The former case is conceptually 
less realistic but attractive since typical computing times per model are a few seconds on modern 
PCs, while NLTE models may require many CPU hours. A third approach, so-called hybrid 
NLTE modelling, can combine the strengths of both types of models, performing sophisticated 
NLTE line- formation computations on LTE atmospheres, at a fraction of the computing cost of 
full NLTE models. This is particularly attractive in those cases where the atmospheric structure 
is close to LTE as for the parameter space investigated here which we will show in the following. 

For our test we have chosen results from the most common among the many available model 
atmosphere codes: on the one hand line-blanketed LTE atmospheres and synthetic spectra 
computed with the Atlas9/Synthe suite of codes [H[2] and on the other line-blanketed NLTE 
models calculated with the Tlusty/Synspec package [3]. Model structures, fluxes and detailed 
synthetic spectra were adopted from published grids, Castelli's ODFNEW models [SJ, the Padova 
spectral library [9], and the OSTAR2002 and BSTAR2006 grids [TO1EEI]. 



1 More luminous stars, either mid/early O-type stars or OB-supergiants close to the Eddington limit, require 
consideration of their stellar winds, i.e. hydrodynamic model atmospheres (and spherical geometry), see e.g. [7J. 



Table 1. Model atoms used with our Detail/ Surface computations. 
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Our own computations rely on Atlas9 atmospheres and NLTE line-formation computations 
with updated and extended versions of Detail/ Surface [HE]. The former solves the coupled 
statistical equilibrium and radiative transfer equations, employing the accelerated lambda 
iteration scheme of [12] . while the latter computes the synthetic spectra, using refined line- 
broadening theories. Model atoms according to Table [T] were employed. More details on our 
hybrid NLTE approach (labeled 'Ads' henceforth) can be found e.g. in [24] . 

4. Model tests 

In order to facilitate the use of existing grids for the model comparisons, a two-step strategy has 
to be followed. In the first step, we compare our tailored Ads modelling with observation 
for strategic points in the investigated parameter space to test its ability to reproduce 
observation. Then, in the second step, Ads models for generic parameters are compared 
to the Tlusty/Synspec and Atlas9/Synthe models at a grid point closest to the real 
star's parameters, thus facilitating a meaningful assessment of the predictive power of all these 
modelling implementations. While the synthetic spectra were convolved with the same rotational 
and macroturbulent broadening profiles as for the comparison with observation, the Gaussian 
('instrumental') profile had to be adjusted to account for the lower resolution of the Padova grid. 
Due to space restrictions we will discuss the case of the well-studied standard star 7 Peg (B2 IV) 
in detail in the following, and comment on more general aspects of the comparisons briefljQ 

For 7 Peg (HD886) we derived T e fr = 22 000K, a surface gravity of log g = 3.95 (cgs units), 
microturbulence £ = 2km/s, solar helium abundance and metal abundances typical for the nor- 
mal early B-type star population in the solar neighbourhood |26t 128] , The nearest parameter 
combination covered by the published grids coincides with these, except that log g = 4.00 and 
solar abundances according to |29j . We have adopted the generic grid parameters nevertheless 
for the Ads model in the second step of the comparison in order to avoid any systematic bias. 

Model atmosphere structures and spectral energy distributions. We begin testing the models 
with a comparison of atmospheric structures in Fig. [2] The temperature structures of the LTE 
and NLTE models agree very well, except for the outermost layers. Differences amount to up 
to 1-2% for all depths relevant for the formation of the spectral lines and continua from the 
far-UV through the optical to the near-IR. Interestingly, the differences are smaller for the low- 
metallicity comparison at l/10th of the solar value, despite the fact that NLTE effects should 
be more pronounced because of increased photon mean-free paths and a harder radiation field. 
A possible explanation of this unexpected finding could be the differences in the metal line 
opacities (which matter more at higher metallicity ). Many more elements are considered in 

2 Our observational sample comprises a total of 29 apparently slowly-rotating B3 to BO dwarf and giant stars, 
for which spectra of high signal-to-noise ratio (S/N ~ 250-800), high resolution (R = AA/A ~40-48 000) and wide 
wavelength coverage were obtained using state-of-the-art Echelle spectrographs on 2m-class telescopes, see [2S] 
for a brief overview. For this sample we have shown that combining a sophisticated analysis technique with our 
Ads modelling allows practically the entire observed spectra in the optical for all these objects 26, 27, 28 to be 
reproduced to similar or even better quality as shown for 7 Peg in the following. 




Figure 2. Comparison of model atmosphere structures and resulting spectral energy 
distributions. Panel (a) presents a comparison of Atlas9 (red) and Tlusty (black) temperature 
structures as a function of mass scale, for models at T e g = 22 000 K, log g = 4.00, £ = 2km/s at 
solar (full lines) and l/10th solar metallicity (dashed lines, shifted by +5000 K). Electron density 
structures are shown in the inset. Panel (b) displays the resulting differences between the Atlas9 
and Tlusty temperature structures (in percent). Panel (c) and (d) compare the Tlusty with 
Atlas9 and Atlas9+Detail Eddington fluxes, respectively. The insets highlight the regions 
around the Lyman and the Balmer jumps. Colour coding according to the legends. 
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Figure 3. Comparison between our computed SED (red line) and observation for 7 Peg (B2 IV): 
IUE spectrophotometry (black line) and UBVRIJHK photometry (squares). 

the Atlas9 than in the Tlusty model, resulting in a stronger line-blanketing effect, exactly as 
indicated in Fig.[2](a). Very good agreement is found also for LTE and NLTE density structures. 
Overall, close agreement of the LTE and NLTE structures is found over the entire parameter 
range considered, for the depths relevant for the formation of the observable line spectra and 
continua, at —2 < log m < —0.3. 

As temperature and density stratifications are not directly observable, we verified the effects 
on line profiles of hydrogen, helium and metal lines by performing line-formation computations 
with Detail & Surface on Atlas9 and Tlusty atmospheres. Indeed, only small differences 
are found for the resulting profiles, which are irrelevant for the stellar parameter and elemental 
abundance determination, see e.g. Fig. 10 of [24] . Maximum effects on metal abundances derived 
from a few individual lines are about 0.05 dex, and less for mean elemental abundances. 

Overall good agreement is also found for the LTE and NLTE spectral energy distributions 
(SEDs), as exemplified in Fig. [2] (c) and (d). The differences in the continua are because of 
the temperature differences at the formation depths. Significant discrepancies occur only in the 
Lyman and helium continua, but it is difficult to test the models in this respect as interstellar 
absorption prevents observations of the extreme- UV for almost all early- type stars. The model 
differences in the extreme-UV vary with T e g. Very good agreement over the entire spectral 
range is achieved at higher T e g, see Fig. 9 of [Mj- Our models can reproduce observed SEDs 
and therefore the global stellar energy output very well, if the observations are corrected for the 
appropriate amount of interstellar reddening. An example, for 7 Peg, is shown in Fig. [3j 

We conclude that LTE and NLTE model atmospheres are essentially equivalent for dwarf and 
giant stars over the range 15 000<T c r< 35 000 K, for most practical applications. 

Line spectra. Resolved line spectra are the most valuable observables as they give information 
about the temperature and density stratifications within the stars' atmospheres. They are also 
well-suited for the determination of stellar parameters and elemental abundances. 

H fc He : Hydrogen and helium provide the deepest features in the spectra of early B-type stars, 
thus sampling the widest extent of the atmospheric stratification. Our overall best-fit Ads 
model is compared to observed profiles for many diagnostic lines of 7 Peg in Fig. [4j We note 
explicitly, that one synthetic spectrum is displayed in all the individual panels - also for all 
following comparisons - and not best fits to the individual lines. Overall, an excellent match of 
observation and theory is obtained. In particular the line wings agree well, being formed in the 
deeper layers of the atmosphere. Small discrepancies are present for the line cores, in particular 
for the stronger He 1 lines, which originate at shallower atmospheric depths. And, apparently 
the treatment of the Hei forbidden components may be improved in some cases. 
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Figure 4. Comparison of strategic Hi and Hei/n in the observed spectrum of 7 Peg (black 
dots) with our global best-fit Ads model (red line). 



About the same fit quality can be achieved throughout the entire parameter space investigated 
here, with additional minor complications arising at the hotter end of the T e g-range, where weak 
stellar winds start to impact Ha and the He 11 A4686 A line, see |24| for examples. 

We conclude that the NLTE modelling with Ads can reproduce the observed Hi and Hei/11 
lines reliably, except for some fine details. Ads models thus provide a good reference for the 
second step of our test, comparison with the Tlusty/Synspec NLTE and Atlas9/Synthe 
LTE models, which is shown in Fig. [5] Both NLTE approaches produce closely resembling profiles 
for this test case, and the good agreement (except for some minor details) holds throughout 
the entire parameter range considered here for the Hi, He 11 and Hei triplet lines. However, at 
higher temperatures the He 1 singlet lines computed with Tlusty/Synspec become significantly 
weaker than the Ads profiles (and the observations) - the well-known singlet-triplet problem [30J . 
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Figure 5. Same as Fig. [4j comparing line profiles computed with Ads, Tlusty/Synspec and 
Atlas9/Synthe (line designations according to the legend). 



Good agreement between the NLTE models and the Atlas9/Synthe profiles is found for the 
Balmer lines and a few sharp Hei lines, for this test case. However, problems with inadequate 
line broadening - Voigt profiles with constant Stark broadening parameter and no account for 
forbidden components - are indicated for the majority of the He I features (in particular all 
diffuse ones) computed with Atlas9/Synthe. This can be resolved simply by implementing 
the appropriate broadening tables. More fundamental are the discrepancies in the depths of 
lines such as HeiAA4921, 5875 and 6678 A, which are genuine NLTE effects. 

The comparison over the entire parameter range reveals severe shortcomings of LTE modelling 
of observed hydrogen and helium lines. Over the T e g-range of 22 000 to 25 000 K, NLTE effects 
quickly set in and affect the cores and wings of the Balmer lines, and subsequently all helium 
lines. Maximum effects are found for highest T e g: the Balmer lines computed in LTE have 
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Figure 6. A): comparison of observations for Nil lines in 7 Peg with our global best-fit Ads 
model. B): comparison of Ads, Tlusty/Synspec and Atlas9/Synthe models for the same 
lines. See figure legends and the text for details. 



equivalent widths about 30% lower than in NLTE, and up to a factor of more than 2 for the 
He 1 lines. Surface gravities determinations based on Balmer-wing fitting is thus subject to 
systematic error, e.g. log g values are overestimated by up to 0.2 dex from LTE analyses of H7. 

Overall, LTE modelling of the Balmer wings and selected Hei lines appears to be largely 
unbiased by systematic error for main sequence stars up to T e ff ~22 000K. Preference should 
then be given to the He I AA3867, 4121, 4437, 4713, 5016 and 5048 A transitions, which are least 
affected by NLTE effects. Either full or hybrid NLTE models have to be used for quantitative 
analyses of stars at higher T e g- when systematic bias is to be avoided. Our Ads approach improves 
on the Tlusty OSTAR2002/BSTAR2006 grids by avoiding the Hei singlet-triplet problem. 




Nitrogen: We discuss here only spectral lines of the N n ion, which in contrast to N in lines are 
observable throughout the entire parameter range. A comparison of important diagnostic lines 
in the optical spectrum of 7 Peg with our global best-fit Ads model is shown in Fig. [6] The fit 
quality is excellent, in total 40 Nil lines were analysed to derive the mean nitrogen abundance. 
Therefore, the Ads model provides a good reference for the comparison with the other codes, 
which is shown in the bottom half of Fig. [6j 

In many cases good agreement is found between the Tlusty/Synspec and Ads NLTE 
line profiles. However, a significant number of differences exists, e. g. for NnAA4630, 5045, 
5679 and 5711 A from the examples shown in Fig. [6j where the Tlusty/Synspec models 
are too weak compared to observation. Abundance analyses based on such features would 
therefore significantly overestimate the nitrogen abundance. Both agreement and discrepancies 
vary throughout the entire parameter range considered here, so that it is difficult to be 
quantify them in brief. In general, the differences increase towards lower gravities. The 
reasons for the differences in the two NLTE approaches result almost entirely from the model 
atom implementations. We tested this by adopting the Nil model atom used in the Tlusty 
computations for Detail/Surface and could tightly reproduce the Tlusty/Synspec solution. 

The comparison of the NLTE Ads reference with the LTE Atlas9/Synthe model shows 
that NLTE effects are often small, except for lines such as NnA4176A from the examples. In 
fact, Atlas9/Synthe profiles are closer to the Ads profiles (and therefore observation) than 
the Tlusty/Synspec profiles in most cases. The NLTE effects become more pronounced with 
decreasing surface gravities. 
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Figure 8. As Fig. [HI for Mgn, Sim/iv and Feli/m lines. 



Oxygen: We investigated the behaviour of O I and O n lines, with similar results to those for 
nitrogen. Because of this we omit a more detailed description here. Noteworthy are the overall 
relatively small NLTE effects, again with the exception of several multiplets, e.g. the near-IR 
triplet OiAA7771-5A. For these lines, large differences between the two NLTE solutions are 
seen, again probably because of the different model atom implementations. 

Carbon : Singly and doubly-ionized carbon were investigated. Again, the Ads model provides a 
good reference for further comparison with other models, as it reproduces observation well, not 
only for 7 Peg (see examples in Fig. [7]) but throughout the entire early B-type main sequence [16] . 

In particular C II was known as one of the most challenging cases for NLTE line formation 
computations for decades. The comparison of Ads with the Tlusty/Synspec solution shows 



disagreement in practically every line. Usually the Tlusty/Synspec profiles are too shallow, 
and in some cases even turn into (unobserved) emission (CnA6462A). The situation is much 
better for C III, except for the strong line C in A4162 A. The problems with C n persist for other 
parameter combinations. Certainly, the Tlusty model atom for carbon needs to be improved 
to allow for meaningful abundance determinations. For the moment, C in should be preferred 
for abundance work based on the BSTAR2006 and OSTAR2002 grids. 

The comparison with the LTE profiles shows that most of the stronger, but also several of the 
weak lines show pronounced NLTE effects, throughout the entire parameter range. An exception 
is CiiA5145A (and its other multiplet members), which stays in LTE for all cases of interest. 
NLTE effects on the Cm lines increase significantly for higher T e g. 

Magnesium, silicon and iron: Examples for the comparison of our global best-fit Ads model 
with observations for important diagnostic Mgn, Sim/iv and Fell/111 lines in 7 Peg are shown 
in Fig. [8} The Mgn and Fell/111 are reproduced reliably. We are at present incorporating a new 
silicon model atom, in particular to improve the fits to several lines not shown here, but the 
old model [21] reproduces lines such as those of the Sim AA4552-4575 A triplet or SiivA4116A 
reasonably well, providing a valuable diagnostic for T e g -determination via the ionization balance. 

Both NLTE solutions, with Ads and with Tlusty/Synspec, agree well for Mgn, but 
disagree for almost all Si and Fe lines. This is of concern in the case of silicon as use of the 
ionization equilibrium would lead to systematically different solutions for T e g- and log g with the 
BSTAR2006 and OSTAR2002 grids, probably inconsistent with other diagnostics. Moreover, 
the Si in triplet is often used to constrain the microturbulence velocity, which therefore may 
also become subject to systematic offsets when based on the Tlusty grids. The iron lines are 
systematically shallower in the Tlusty/Synspec models, implying the derivation of higher iron 
abundances in applications. Because of the role of iron as the main line-opacity source this has 
further effects on atmospheric structure computations and therefore also on the atmospheric 
parameter determination. It is highly important to improve these model atoms for Tlusty. 

While the LTE and NLTE solutions for Mgll are very similar for the parameters of 7 Peg 
(Mg 11 A4390 A apparently differs because of its location in the wing of an ill-fitted He 1 line) , the 
iron lines are weakly and most of the silicon lines strongly NLTE affected. LTE modelling of 
silicon lines in early B-type stars can lead to systematic errors in atmospheric parameters when 
using ionization equilibria, and to unreliable elemental abundances. LTE iron abundances are 
in general somewhat underestimated. 

5. Conclusions and Recommendations 

We have seen that the data from the Gaia mission will need to be complemented by ground- 
based (optical) spectroscopy for early-type stars in order to put meaningful constraints on 
atmospheric parameters and elemental abundances. Three approaches for the modelling of stellar 
atmospheres in the T e g-range of 15 000 to 35 000 K, corresponding to spectral types of about B4 
to 09, have been tested here: full LTE (with Atlas9/Synthe), full NLTE (Tlusty/Synspec) 
and hybrid NLTE (Ads). The aim of the work was to determine the best means, presently 
available, to analyse optical spectra of such early-type stars, which is of course of much broader 
relevance than for Gaia-related science alone. 

The comparison showed that pure LTE modelling can yield meaningful results for the low- 
T c fj regime, up to about 22 000 K on the main sequence, if the right spectroscopic indicators 
are employed. An overall good global match between synthetic and observed spectra cannot be 
expected to be obtained as many lines show NLTE effects. 

Pure NLTE modelling with the available Tlusty/Synspec grids is conceptually superior. 
However, NLTE atmospheric structures and SEDs turn out to be close to LTE in the parameter 
space considered. NLTE line spectra of ions with 'simple' atomic structures, like H 1, He 1/11 (with 
some restrictions at higher T e g) or Mgll turn out to be reliable. However, oversimplified model 



atoms cause problems for many of the more complex ions/elements. Blind use of the synthetic 
spectra published in the OSTAR2002 and BSTAR2006 grids should therefore be avoided. The 
fact is that NLTE computations based on inappropriate model atoms can potentially introduce 
larger systematic errors to quantitative analyses than simple LTE modelling. The Tlusty model 
atoms certainly need extensive revision before a good global match with observations can be 
obtained with this package. 

Our NLTE line- formation computations with Detail/ Surface on Atlas9 atmospheres 
combine the advantage of low to moderate CPU time requirements with realistic model atoms, 
facilitating a robust global match of observed spectra to be obtained, see also [271 121] • It 
is reassuring that the same model atoms provide as good a match of the synthetic spectra 
to observations also for other kinds of stars, with drastically different chemical composition 
such as subdwarf OB stars [31] or extreme helium stars (32J, [331 EH- First applications of our 
models combined with a sophisticated analysis methodology promise breakthroughs for several 
astrophysical fields, see e.g. [25] . 
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